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Abstract 

We discuss the elliptic flow dependence on pseudorapidity and number of participating 
nucleons in the framework of string percolation, and argue that the geometry of the 
initial overlap region of interaction, projected in the impact parameter plane, determines 
the experimentally measured azimuthal asymmetries. We found good agreement with 
data. 



^^' The discovery of the large elliptic flow V2 was one of the most important achievements 

D , at RHIC experiments [1-9]. A non vanishing anisotropic flow exist only if the particles 

measured in the final state depend not only on the local physical conditions realized at the 
production but as well on the global event geometry. In a relativistic local theory, this non 
local information can only emerge as a collective effect, requiring interactions between 
^ \ the relevant degrees of freedom, localized at different points of the collision region. In 

this sense, anisotropic flow is particularly unambiguous and convincing manifestation of 
p^ , collective dynamics in heavy ion collisions [10]. 

l/^ ' The elliptic flow V2 can be qualitatively explained as follows. In a collision at high en- 

ergy the spectators are fastly moving opening the way, leaving behind at mid- rapidity an 
f— ^ ' almond shaped azimuthally asymmetric region of QCD matter. This spatial asymmetry 

f^ , implies unequal pressure gradients in the transverse plane, with a larger density gradient 

in the reaction plane ( in- plane). As a consequence of subsequent multiple interactions 
between degrees of freedom this spatial asymmetry leads to an anisotropy in momentum 
plane. The final particle transverse momentum is more likely to be in- plane than in the 
. , out- plane, with ^2 > as predicted [11]. 

j^ ' The basic idea of our model [12] is that the angular azimuthal anisotropy associated 

to the geometry of the first stages in the collision - the projected almond- influences in a 
determinant way the presence or not of the flow. In other words if the projected overlap 
region was a circle we would have W2 = 0. As in the almond case the small axis is in the 
reaction plane, corresponding to higher matter density, then V2 > 0. 

Our model was introduced in [12] and a discussion of applications and conjectures 
were presented dependence of V2 on the produced hadron, validity of quark counting 
rules, applications to nuclear reduction factors, etc. Here we just want to call attention 
to V2 as a function of pseudorapidity rj and the number of participants Na for nucleus A, 
after the integration over p^. The triangle shape shown by the data on the dependence 
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of V2 on r] it is not easy reproduced by models as it has been recently emphasized [13]. 
We show that string percolation model is able to do it. 

The string percolation model [14] develops around the concept of transverse density 



l' - C-^)'N'^ (1) 

, where N'' is the number of longitudinal strings formed in the collision, vq is the radius 
of the single string and R the effective radius of the interaction overlap region S in the 
impact parameter b, 

S = 7ri?2 (2) 

with 



S = 2R\[cos-^iP)-^^l-p^] (3) 

, Ra being the nuclear radius and 

Two relations, one for the particle density dn/drj and the other for the average trans- 
verse momentum squared define the essential features of the model [13,14] 

dn/dri = F{rj*)N'')j,i (5) 

and 

<p''t>=<Pt>i/F{v') (6) 

where /ii and < p^ >i are single string parameters and F[rf) is the colour reduction 
factor [15]: 

/l - e-')' 

We introduce now two reasonable approximations: that TV* is proportional to the number 
of binary collisions and that R is proportional to the proton radius, 

N' ~ N^N^f (8) 

and 

R = RpN\^\ (9) 

where N^ and Rp are proton parameters and Na is the number of participants from 
nucleus A. 

From (1), (8) and (9) we obtain 

V%. = vlN'f. (10) 

By using (9) and (1) on (5) one obtains 

1 d^ f(^*)^*(^)2^^ (n) 



jV2/3 (Ijj ro 

and we observe that the right hand side of (11) and (6) are deeply related. 
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This kind of results appears in the Color Glass Condensate (CGC) [16] and in string 
percolation [14,17]. Note that (11) can be written in the form 

Trrn , 1 dn. 



v/(l-exp-ry*V = — I o^tI (^2) 

This relation, as we shall see, is essential to understand the (pseudo)rapidity and 
number of participants per nucleus dependence of W2: V2{r],NA)- Note that in (12) small 
rj* corresponds to large 77 and large 77* to small 77. 

Regarding p^ distributions, we started with Schwinger gaussian formula, including 
fusion and percolation (via F{rf)) and clustering fluctuations (via the parameter k{rf)) 
to obtain[18]: 

(Pn _dnk-l F{r]*) 1 



dp'^drj dr] k < _p|. >i n 1 Zil^Mi^fe ' 



(13) 



fe<P5,>i • 

Most of the RHIC data are well described by formula (13) [12,18,19]. 

In order to discuss directional production along the azimuthal angle 1^, we shall 
introduce a convenient variable 

X - F{7^')pI, (14) 

and Xip 

X^ = F{Tf^)pl. (15) 



with 






Vt^v\ — f (16) 



such that we can simplify notation 

dn 



dp\dy 



dn 



dp^dydcp 
Expanding now X^ or i?^ around X or R^ we write 



fix) (17) 



nx^) (18) 



^^'--■fi^)i dR 



fiX^) ^ z/(^)[ ;o^^ (^^ - R')] (19) 



Note that (12) satisfies the normalization condition 

f{X^)d^ = f{X) (20) 

because R^ —< R^ > [12]. Finally we obtain for V2, a function of several variables 
including p^, rj and Na, 

2 f/' ,„ ,,fl/,l e-''-f(ry*)^ F(r;*)p|/ < p^ >i 
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which we shall write as the product of tree factors, 

V2 ^ MnF{rj')p'^], (22) 



2 /""/^ R ^ 

- / d^cos{2^){-^ ), (23) 

TI" Jo -K 



or, having present that 



R^ sin ip — a 



(24) 



and 



Ra sin (p 

a — sin" {f3sin(p), (25) 

M = - / ^' rf^££f(Mf!!!!(^^^(:^)2, (26) 



,ti _ 1 e** - F(r/*)2 
2 F(r/*^2 



^1 - r ..: r (27) 



and 



' fe<p^>i 

Let us next look at the factor [(/j], (23) in (21) and consider, for fixed rj and y/s, two 
limits: 

i) 5 ^ or /? ^ or iV^i ^ A which implies, (25), a ^ 0. We then have [ip] — > 0, 
[r/*], (27),^ constant, and [F(r/*)p|,] -^ 0, or: 

7Va^A,«2(Pt)^0 (29) 



a) b -^ 2Ra or /3 -^ 1 or Na — > which implies, (18), a — > (p. We then have [ip] -^ 0, 
[77*] — > constant, [F{ri'^)p^] — > some finite function of p^, or: 

Na^O, V2{Pt) ^ (30) 

If wc look now to the p^ dependence of V2 in [F{ri*)p'^], (28) we see that U2 — > as 
p^ ^ and V2 ^ k ^ constant, as p^ —> 00. This is observed in data [20]. 
We perform next the integration in p^, weighted by ^ "^"^ / ^ , to obtain: 

, [?7*]' being different from [rj*] and 77* being related to rj by relation (12). Applying now 
to V2{r]) the arguments used for V2{px) we have: 

i)Na -^ A, V2{v) -^ 0, (32) 

with [?7*]' being some negative number depending on 77; 

ii)NA -^ 0, V2{v) -^ 0, (33) 
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with [77*]' -> 0. 

In conclusion, both V2{px) and V2{ri) go to zero as Na — > A and Na -^ 0. In the case 
of W2(p^) with p^ = 0, V2{p'^) is identicahy zero. 

As [ry*] is, in modulus, a growing function of 77*, it is clear that V2{ri), at fixed 77, is a 
growing function of energy and of Na see [20] . 

Regarding V2 normalized by the eccentricity e, 



, with 

having the limits 
and 



V2{ri)h. (34) 

P^Q,e^ P^O (36) 

/? ^ 1 , e ^ 1 (37) 



We see that 

i) Na ^ ^, i? ^ 0, V2(rj)/e — > constant, increasing with y^ and Na 

a) Na-^ 0,13^1, V2{v)/e -^ 0. 

In order to compare with experimental data the dependence of V2 on the pseudora- 
pidity, we start with the ^ data of PHOBOS collaboration [20] taken at Npart = 211. 
From formula (12) we compute 77* at each value of 77 and then V2 using equation (31). Our 
result together with the experimental data [20] is presented in fig 1. In the same way, 
using equation (31) we compute the dependence of V2 on the number of participants. In 
fig 2. we show our results together with the experimental data. In both cases, rapidity 
and centrality dependence, the agreement is very good. 

Summarizing up, the analytical formulae (21) and (34) obtained in the framework 
of string percolation are able to describe rightly the dependence of the elliptical flow on 
rapidity and centrality. 

1. Acknowledgements 

This work is under the project FPA2008-01177 of Spain, and of the Xunta de Galicia. 
J. Dias de Deus thanks the support of the FCT/Portugal project PPCDT/FIS/575682004. 

References 

[1] K. Adcox et al. [PHENIX Collaboration], Nucl. Phys. A 757 (2005) 184. 

[2] J. Adams et al. [STAR Collaboration], Nucl. Phys. A 757, 102 (2005). C. Adlcr et al. [STAR 

Collaboration], Phys. Rev. Lett. 87, 182301 (2001). 
[3] S. Manly et al. [PHOBOS Collaboration], Nucl. Phys. A 774, 523 (2006). 
[4] N. Borghini and U. A. Wiedemann, J. Phys. G 35, 023001 (2008). 
[5] J. Y. OUitrault, Phys. Rev. D 46, 229 (1992). 
[6] P. Huovinen, P. F. Kolb, U. W. Heinz, P. V. Ruuskanen and S. A. Voloshin, Phys. Lett. B 503, 58 

(2001). 
[7] L. Bravina, K. Tywoniuk, E. Zabrodin, G. Burau, J. Bleibel, C. Fuchs and A. Faessler, Phys. Lett. 

B 631, 109 (2005). 

5 





- 
























0.05 
0.04 
0.03 


■ 


+ 


■+ 


+ 


■+ 


+ 


■+ 


+ 


+ 


+ 


+ 




0.02 


L 






















I 


0.01 


'- 








1 










1 







-4 -3 -2 



-1 ^0 1 



Figure 1: Elliptic flow as a function of pseudorapidity for Npart = 211 in Au+Au collisions at energy 
^/s = 200 GeV. Dots in blue are used for our results and bars in red are data taken from reference [20]. 




Figure 2: Elliptic flow dependence on the number of participants, at energy ^/s = 200 GeV. Results 
compared to PHOBOS data. Linos in blue arc used for our results and red lines are data taken form 
reference [20]. 
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